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Determining the global minimum of Higgs potentials via Groebner bases 

- applied to the NMSSM 
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Determining the global minimum of Higgs potentials with several Higgs fields like the next-to- 
minimal supersymmetric extension of the Standard Model (NMSSM) is a non-trivial task already 
at the tree level. The global minimum of a Higgs potential can be found from the set of all its 
stationary points defined by a multivariate polynomial system of equations. We introduce here the 
algebraic Groebner basis approach to solve this system of equations. We apply the method to the 
NMSSM with CP conserving as well as CP violating parameters. The results reveal an interesting 
stationary-point structure of the potential. Requiring the global minimum to give the electroweak 
symmetry breaking observed in Nature excludes large parts of the parameter space. 



1. INTRODUCTION 

It is a non-trivial task to find the global minimum for 
Higgs potentials with a large number of Higgs fields. For 
instance in the next-to-minimal supersymmetric exten- 
sion of the Standard Model (NMSSM) Q, the Higgs sec- 
tor consists of two (complex) electroweak doublets and 
one (complex) electroweak singlet, that is 8 real fields 
from the doublets plus 2 real fields from the singlet. The 
conventional approach based on the unitary gauge re- 
quires the global minimum to be found in a 7-dimensional 
field space. 

In this paper we introduce an algebraic approach to 
determine the global minimum of the Higgs potential. 
We describe how to compute all stationary points, from 
which then the one with the lowest value of the poten- 
tial can be identified as the global minimum. We apply 
the method to the NMSSM, where we can reveal a quite 
surprising structure of stationary points, that is minima, 
maxima, and saddle points with different behaviour with 
respect to the symmetry breaking of the SU(2)l x U(1)y 
electroweak gauge group. 

The global minimum of the Higgs potential gives the 
expectation values of the Higgs fields at the stable vac- 
uum. Parameter values for the Higgs potential are thus 
considered acceptable only, if the global minimum of the 
Higgs potential occurs for Higgs field vacuum expecta- 
tion values, which induce the spontaneous breakdown of 
SU(2)l x U(l)y to the electromagnetic U(l) em at the 
observed electroweak scale v w 246 Ge V. 

We consider the tree level Higgs potential for general 
models with two Higgs doublets and an arbitrary num- 
ber of additional Higgs singlets. The first step is to no- 
tice that the potential is restricted by renormalisability 
and gauge invariance. Renormalisability requires at most 
quartic terms in the real fields in the potential. Elec- 
troweak gauge invariance restricts the possible doublet 
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terms in the potential, since only gauge invariant scalar 
products of doublets can occur. Substituting the dou- 
blet fields by appropriate functions of these invariant 
terms, we eliminate all gauge degrees of freedom from 
the potential and effectively reduce the occurring powers 
in the doublet terms. The method to base the analysis 
on quadratic gauge invariant functions was introduced 
already in context with the general two-Higgs doublet 
model HQ. 

If the potential is bounded from below, the global min- 
ima are given by the stationary points with the lowest 
value of the potential. The stationarity conditions form a 
non-linear, multivariate, inhomogeneous polynomial sys- 
tem of equations of third order. In this work we want 
to introduce a systematic approach to solve these - in 
general quite involved systems of polynomial equations. 
Since we eliminate the gauge degrees of freedom by means 
of gauge invariant functions the number of complex solu- 
tions is assumed to be finite, which implies the absence 
of further continuous symmetries. We propose to deter- 
mine the stationary points by a Groebner basis compu- 
tation, which is well established in ideal theory 0, H, l1 • 
The Groebner basis was originally introduced to solve 
the ideal membership problem. Constructing this Groeb- 
ner basis in an appropriate order of the monomials (the 
terms of the polynomials including coefficients), for in- 
stance the lexicographical ordering, and subsequent tri- 
angularisation allows to solve the initial system of equa- 
tions. The method guarantees that all stationary points 
are found. 

We apply the method introduced to the NMSSM. For 
the computation of Groebner bases as well as the sub- 
sequent steps to solve the systems of equations we em- 
ploy the freely available open-source algebra program 
SINGULAR |2|. We find that large parts of the param- 
eter space of the NMSSM Higgs potential can be ex- 
cluded by requiring the global minimum to lead to the 
electroweak symmetry breaking observed in Nature. We 
illustrate this by determining the allowed and forbidden 
ranges for some generic parameter values for this model. 

In the literature we found many conditions which con- 
strain parts of the parameter space in the NMSSM, see 
for instance 0, 0, . Typically the conditions used are 



2 



necessary but it is not clear if they are also sufficient 
to ensure that the resulting theory is acceptable. We 
also want to mention that there is a purely numerical 
approach to determine the global minimum in the effec- 
tive one-loop NMSSM Higgs potential The aim of 
our present work is to systematically reveal all stationary 
points of the tree level potential by solving the system of 
equations which originates from the stationarity condi- 
tion. We show, that this can be done for the full param- 
eter space including CP violation. With our method we 
can decide unequivocally if a given parameter set of the 
tree level Higgs potential leads to an acceptable theory 
or not. 



the doublet degrees of freedom enter solely via products 
of the following form: 



'Pi'Pj 



with i, j <G {1, 2}. 



(2.5) 



It is convenient to discuss the properties of the potential 
such as its stability and its stationary points in terms 
of these gauge invariant quadratic expressions. This was 
discussed in detail for THDMs and also extended for the 
case of more than two doublets in . We recall the main 
steps here. 

We arrange all possible SU(2)l x U(1)y invariant 
scalar products into the hermitian 2x2 matrix 
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K ._ Mvi <Pa<Pi\ 



(2.6) 



We consider the tree-level Higgs potential of models 
having SU(2)l x U(1)y (weak isospin and hyperchargc) 
electroweak gauge symmetry. In particular we study 
models with two Higgs doublets and n additional real 
Higgs isospin and hyperchargc singlets. This includes in 
particular THDMs, where we have no additional Higgs 
singlets, and the NMSSM with one additional complex 
Higgs singlet corresponding to two real singlets. We as- 
sume both doublets to carry hypercharge y = +1/2 and 
denote the complex doublet fields by 



(Pi(x) 



1,2. 



(2.1) 



For the singlets we assume real fields which we denote by 



1, 



(2.2) 



We remark that in supcrsymmctric models like the 
NMSSM the two Higgs doublets Hd, H u carry hyper- 
charges y — —1/2 and y = +1/2 respectively. This can 
be translated to the convention used here by setting 



r-2 



where 
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(2.3) 



(2.4) 



and consider its decomposition 
1 



A"; 



A'q^ 



K a cr?-, 



(2.7) 



where a a are the Pauli matrices. The four real coefficients 
in this decomposition are 

K = d<Pi, K a = a =1,2, 3, (2.8) 

where here and in the following summation over repeated 
indices is understood. The matrix 1|2.6|) is positive semi- 
definite, which implies 

A > 0, Aq - K\ - Af - Af > 0. (2.9) 

On the other hand, for every hermitian 2x2 matrix A„- ■ 
satisfying 1)2. 9|1 there exist fields l|2.1(l satisfying 1)2. 6[l . 
see [3|. It was also shown in Q that the four quantities 
Ko,K a satisfying (|2.9|) parametrise the gauge orbits of 
the Higgs doublets. Using the inversion of (12.811 . 



(pi?! = (A + A 3 )/2, ip\<p 2 = (Ai + iK 2 )/2, 
<pl<(>2 = (A - A 3 )/2, <4<^ = (Ai - iK 2 )/2, 



(2.10) 



we can replace the doublet terms of the potential - due 
to renormalisability - by at most quadratic terms in the 
real functions Ao, Ai, K 2 , and A3, which simplifies the 
potential and eliminates all SU(2)l x U(1)y gauge de- 
grees of freedom. We thus end up with a potential of the 
form V(K , Ai, A 2 , A' 3 , <f> lt . . . , 



Complex singlet fields are embedded in our notation by 
treating the real and imaginary parts of the complex sin- 
glets as two real singlet degrees of freedom. 



2.1. Gauge invariant functions 

In the most general SU(2)l x U(1)y gauge invariant 
Higgs potential with the field content described above, 



2.2. Stationary points 

To determine the stationary points of the Higgs poten- 
tial we consider V(Kq, Ki,K 2 , A3, <f)\, . . . , 4> n ) and take 
the constraint l|2.9|l into account. We distinguish the pos- 
sible cases of stationary points by the SU(2)l X J7(l)y 
symmetry breaking behaviour which a vacuum of this 
type would have (see Q): 
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• unbroken SU{2) L x U(1) Y 
with 



A stationary point 



K =K 1 =K 2 = K 3 = 0. 



(2.11) 



A global minimum of this type implies vanish- 
ing vacuum expectation value for the doublet 
fields lIU and therefore the trivial behaviour with 
respect to the gauge group. The stationary points 
of this type are found by setting all Higgs-doublet 
fields (or correspondingly the Kq as well as the K a 
fields) in the potential to zero and requiring a van- 
ishing gradient with respect to the remaining real 
fields: 



v v(4>i 



= o. 



(2.12) 



• fully broken SU{2)l x U(1)y ■ A stationary 
point with 



K > 0, 
Kl - K\ - K% - K\ > 0. 



(2.13) 



A global minimum of this type has non-vanishing 
vacuum expectation values for the charged compo- 
nents of the doublets fields in 1)2. lfl . thus leads to 
a fully broken SU(2) L x U(1) Y , see 0. The sta- 
tionary points of this type are found by requiring a 
vanishing gradient with respect to all singlet fields 
and all gauge invariant functions: 

V V(K , Kx, K 2 , K 3 , </>!,..., <j> n ) = 0. (2.14) 

The constraints i|2.13[) on the gauge invariant func- 
tions must be checked explicitly for the real solu- 
tions found. 

• partially broken SU{2)l x U(1)y '■ A stationary 
point with 



K o >0, 
Kl - Kl - Kl - Kl = 0. 



(2.15) 



For a global minimum of this type follows the 
desired partial breaking of SU{2)l x U(1)y down 
to U(l) em , see @- Using the Lagrange method, 
these stationary points are given by the real solu- 
tions of the system of equations 

V[V{K ,K x ,K 2 ,Kz,ct>x,...,<t>n) 

-u ■ {Kl - Kl - Kl - Kl) } = 0, (2.16) 
Kl - Kl - Kl - Kl = 0, 

where u is a Lagrange multiplier. The inequality 
in (|2.15|) must be checked explicitly for the solu- 
tions found for 12.16fl . 

For a potential which is bounded from below, the global 
minima will be among these stationary points. Solving 



the systems of equations (I2.12|) . (|2.14[) . and (|2.16|l . and 

inserting the solutions in the potential, we can therefore 
identify the global minima as those solutions which have 
the lowest value of the potential. Note that in general 
there can be more than one global minimum point. 

From the mathematical point of view we have 
with l|2.12l 12.141 12.1tjfl to solve non-linear, multivariate, 
inhomogeneous systems of polynomial equations of third 
order. We want to demonstrate that this is possible, 
even if the number of fields is large, like in the NMSSM. 
The most involved case is given by (I2.16|) . which for the 
NMSSM consists of seven equations in seven indetermi- 
nates, namely six real fields and one Lagrange multiplier. 
We stress that at this point we assume the number of 
complex solutions of (|iTT2l ICTH YTTtyi to be finite, that 
is, we have systems of equations with a finite number of 
solutions. The gauge invariant functions avoid spurious 
continuous sets of complex solutions, which we found to 
arise in the case of the MSSM as well as the NMSSM 
if the stability conditions are formulated with respect to 
the Higgs fields (|2.1|l in a unitary gauge. 

The solution of multivariate polynomial systems of 
equations is the subject of polynomial ideal theory and 
can be obtained algorithmically in the Groebner basis 
approach 0. See appendix 1X1 for a brief introduction to 
this subject. Within this approach the system of equa- 
tions is transformed into a unique standard form with 
respect to a specified underlying ordering of the poly- 
nomial summands {monomials). This unique standard 
form of the system of equations is given by the corre- 
sponding reduced Groebner basis. If the underlying order 
is the lexicographical ordering, the unique standard form 
consists of equations with a partial separation in the in- 
determinates. We use a variant of the F& algorithm ^2] 
to compute the Groebner bases. A Groebner basis com- 
putation is generally much faster if the standard form is 
computed with respect to total degree orderings and then 
transformed into a lexicographical Groebner basis. The 
transformation of bases from total degree to lexicograph- 
ical ordering is done with the help of the FGLM algo- 
rithm 0. Finally, the system of equations represented 
by the lexicographical Groebner basis has to be triangu- 
larised. The decomposition of the system of equations 
with a finite number of solutions into triangular sets is 
performed with the algorithm introduced in |l4j . Each 
triangular system consists of one univariate equation, one 
equation in 2 indcterminates, one equation in 3 indeter- 
minates and so forth. This means that the solutions are 
found by subsequently solving just univariate equations 
by inserting the solutions from the previous steps. 

The construction of the Groebner basis as well as the 
triangularisation are done algebraically, so no approxi- 
mations are needed. However, the triangular system of 
equations contains in general polynomials of high order, 
where the zeros cannot be obtained algebraically. Here 
numerical methods arc needed to find the in general com- 
plex roots of the univariate polynomials. 
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In more involved potentials, like the NMSSM, the algo- 
rithmic solution is considerably simplified (or even made 
accessible), if the coefficients of the polynomials are given 
in form of rational numbers. Since rational numbers are 
arbitrarily close to real numbers and moreover the physi- 
cal parameters are given only with a certain precision this 
does not limit the general applicability of the method in 
practice. 

All algorithms for the computation of the Grocbner 
basis with respect to a given order of the monomials, the 
change of the underlying order, the triangularisation, and 
the solution of the triangular systems are implemented 
in the SINGULAR program package Q. The solutions 
obtained can be easily checked by inserting them into 
the initial system of equations. Moreover, the number of 
complex solutions, that is the multiplicity of the system, 
is known, so we can easily check that no stationary point 
is missing. 



3. STATIONARY POINTS IN THE NMSSM 

Now, we want to apply the methods introduced in sec- 
tion □ to the NMSSM. 



3.1. The NMSSM Higgs potential 

The NMSSM Higgs sector contains two doublets and 
one singlet, 

H u = , H d = , S, (3.1) 

with the tree level Higgs potential 
Vnmssm = V f + V d + V soft HE3, where 

V F =\XS\ 2 (\H U \ 2 + \H d \ 2 ) + \XH u H d + K S 2 \ 2 , 
V^J-f{\H d \ 2 -\Hu?) 2 + \ 9 2 \HtH d \\ 

(3 2) 

V soi t =m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 |S| 2 
+ [XA x SH u H d + ^A K S 3 + h.c.]. 

We use the notation H u H d = e a /3(-ff u ) a (i?d)^ = 
H+H- - H a u H* and g = s/g 2 +g' 2 , where g and g' 
are the SU(2)l and U(l)y gauge couplings, respectively. 
The parameters of the potential are given by the experi- 
mentally fixed electroweak gauge couplings and 

^,K,m 2 Hu ,rn 2 Hd ,m 2 s ,A x ,A K . (3.3) 

The quartic terms of the potential l|3.2(l are positive 
for any non-trivial field configuration, if both A, n are 
non- vanishing. The potential is therefore bounded from 
below for all cases considered here, and stability has not 
to be checked any further. 



We translate the NMSSM Higgs potential to the for- 
malism described in the previous section, where we de- 
compose the complex singlet field into two real fields ac- 
cording to S = S re + iSi m . In this notation the potential 
is given by 

V F = J|A| 2 (A 2 + K\ + 4A (5 2 e + S 2 J) 
+ \k\ 2 (S 2 e + S 2 m ) 2 

- Re(A^*) (Ax(5 2 e - S 2 m ) + 2K 2 S re S im ) 
+ Im(A K *) (K 2 {S 2 e - S 2 n ) - 2K x S n S im ) , 

Vd = \fK 2 + i 3 2 (A 2 - A 2 - A 2 - A 2 ) , 

Vsoft = ^m 2 Hu (A - A 3 ) + \m 2 Hd (A + A 3 ) 

+ m 2 (S 2 e + S 2 m ) 
-Re(XA x )(K 1 S re -K 2 S mi ) 
+ lm(XA\) (K 2 S re + KxS im ) 

2 3 2 

+ g R.e(ftA K ) (S re — 3S re S im ) 
+ hm( K A K )(Sf m -3S 2 e S. im ). 

(3.4) 

For given values of the potential parameters 1)3.3(1 we can 
find all stationary points of the NMSSM by solving the 
systems of equations (|2.12ll2~T4ll2.16|l as described above. 

3.2. Choice of parameters 

In order to fix experimentally known parameters like 
the electroweak scale it is inappropriate to choose numer- 
ical values for the parameter set 1(3. 3J) . Instead, we ex- 
press different original parameters in terms of the desired 
vacuum expectation values of the neutral components of 
the Higgs doublets and the Higgs singlet, the mass of the 
charged Higgs boson and a CP -violating phase. To this 
end we introduce the vacuum expectation values of the 
Higgs fields, 

(H d ) = , (3-5) 

(S) = \e**v 8 , 

which are parametrised by v Ul v d , vs, and the phases <p u , 
and (ps- The stationarity condition of the potential at the 
vacuum values for the fields requires a vanishing gradient 
with respect to the Higgs fields. These equations relate 
the soft-breaking mass parameters with the vacuum ex- 
pectation values. As usual, we define v 2 = v 2 ^ + v 2 and 
tan/3 = v u /v d . Further, if we write the complex param- 
eters A, k, A\, and A K in polar coordinates with phases 
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<^A, i5k, <$a a , $a k and introduce the abbreviations 

Sedm = 8\ + Vu + Vs, 5' K = 6 K + 3(fs, (3.6) 

the initial parameters of the potential 1(3. 3|1 can be re- 
placed by the new set of parameters 

A, re, | A K | , tan/3, v s , m H ± , sign R K , <5 E dm , S' K . (3.7) 

plus the electroweak scale v « 246 GeV. Note, that it is 
not sufficient to supply the length |^4 K |, in addition we 
have to fix the sign of 

R K = l=R c ( K A K e^ s ). (3.8) 

In the mass matrix of the Higgs scalars, the CP vio- 
lating entries which mix the "scalar" with the "pscu- 
doscalar" fields are proportional to the imaginary part 
of exp[i(6 EDM - <%)]■ 



3.3. Numerical results 

As a numerical example we choose the parameter val- 
ues 

A = 0.4, re = 0.3, \A K \ = 200 GeV, 
tan/3 = 3, vs = 3v, m H ± = 2v, (3-9) 
signi? K = -, Sedm = 0, 5' K = 

and consider the variation of one parameter at a time 
with the values of the other parameters in 1)3.9(1 kept 
fixed. For a given point in parameter space we com- 
pute all stationary points of the NMSSM potential as 
described above. As mentioned in section the Groeb- 
ner basis construction is performed with numerical co- 
efficients. Here we use a precision of 12 digits for the 
input parameters 1(3.3(1 . which are determined from the 
values for the parameters 1(3.7(1 . The roots of the univari- 
ate polynomials are found numerically, where we choose 
a precision of 100 digits. Our approach allows to use ar- 
bitrary precisions in both cases. We verify that the errors 
of the approximate statements described in the following 
are under control. For generic values of the parameters 
we find 52 complex solutions: 7 corresponding to the un- 
broken, 38 to the partially broken, and 7 to the fully 
broken cases. The number of real and therefore relevant 
solutions depends on the specific values of the parame- 
ters. 

As expected from the Z% symmetry of the potential, 
we find either 1 or 3 solutions sharing the same value of 
the potential within the accuracy of the numerical roots. 
From the computed stationary points only those may be 
accepted as global minima which correspond to the initial 
vacuum expectation values (up to the complex phases), 



that is which fulfil 

(3.10) 

Since for non-vanishing A, re the potential is bounded 
from below, the stationary point with the lowest value 
of the potential is the global minimum. Further, we de- 
termine for every stationary solution, whether it is a local 
maximum, local minimum or a saddle point. For the reg- 
ular solutions, i.e. the partially- and fully-broken cases, 
this is achieved via the bordered Hessian matrix (see 
for instance [3) in terms of K ,K 1: K 2 ,K 3 

, S re , Sim . 

This takes all powers of the doublet fields into account, 
which allows for a definite decision on the type of the 
stationary point also for frequently encountered par- 
tial breaking solutions where at least one mass squared 
is zero and the others have the same sign. For ir- 
regular solutions, i.e. the non-breaking solutions with 
Kq = 0, the Lagrange formalism can not be used since 
the gradients of the two boundary conditions with re- 
spect to Kq, Ki, K 2 , -K3, S re , Sim become linearly de- 
pendent. Instead we resubstitute the original fields 
H+ , H°, H®, H~[ , S re , Si, n in this case and consider the 
free Hessian matrix with respect to these fields. This 
turns out to be sufficient in practice to judge on the type 
of the stationary points. 

Figures and [21 show the values of the potential 
at all stationary points for the parameter values l|3.9(l 
and the cases where successively one of the parameters 
X,K,vs,rn H ±,S' K is varied. Each curve in the Figures 
represents 1- or 3-fold degenerate stationary potential 
values, where the gauge symmetry breaking behaviour 
of the solutions is denoted by different line styles. Ex- 
cluded parameter regions, where the global minimum 
does not exhibit the required expectation values ((3.10|) 
are drawn shaded. As is illustrated by the figures we 
find that substantial regions of the NMSSM parameter 
space are excluded. For some excluded parameter re- 
gions, the partially breaking solutions with the required 
vacuum expectation values Q3.10|) arc saddle points (see 
for instance Fig. ^ top). This means they can be 
discarded as global minima without calculation of the 
other stationary points. However, this is not always 
the case. Obviously from Fig. [21 top, we find an up- 
per bound for vs- For the plotted vs larger than this 
upper exclusion bound the solutions fulfilling 1(3.101) arc 
still pronounced minima, i.e. the mass matrices have 
positive eigenvalues, but they are no longer the global 
minima. The influence of the CP violating phase 5' K 
is shown in Fig. [21 bottom. Note that S' K — > —6' K is 
not a symmetry of the potential. However, the poten- 
tial is invariant under (S' K ,K 2 ,Si m ) — > —(S' K , K 2 , S m ), 
that is (S' K ,H u ,H d ,S) (S' K ,H*,H^,S*), if the resid- 
ual phases are chosen as in (|3.9fl . Therefore the station- 
ary values of the potential in Fig. [21 depend only on \S' K \. 
Also we want to note that in all figures shown here there 
is a non-breaking stationary point with potential values 
slightly above that of the wanted global minimum. We 
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FIG. 1: Values of the NMSSM potential at its stationary points in dependence of k, A. The following parameters are kept 
constant unless explicitly varied: A = 0.4, k = 0.3, \A K \ — 200 GeV, tan /3 = 3, us = 3v,m H ± — 2u,signi? K = — , <5edm = 5' K = 0. 
Each line corresponds to 1 or 3 stationary points sharing the same value of the potential. The different line styles denote saddle 
points, maxima, and minima. The labels 'none', 'full', and 'partial' denote solutions of the classes with unbroken 12.1211 . fully 
broken 12.141 . and partially broken 12.1611 SU(2)l X U(1)y , respectively. For solutions of the partially broken class, it is also 
denoted whether they correspond to the 'required VEVs' v u ,Vd, vs or to 'other VEVs'. Excluded parameter values, where the 
global minimum does not exhibit the required vacuum expectation values, are drawn shaded. 



find that this effect is not coincidental for the parame- 
ters <|3.9[) chosen here, but rather a generic feature of the 
NMSSM. Within the CP conserving parameter range 

Ae]0,l], kg]0,1], A K € i]0,2500] GeV, 

tan/3 e]0,50], v s <g]0,5000] GeV, (3.11) 
m H ± e]0,2500] GeV 

we find a typical separation of the potential values for 
these non-breaking solutions and the solutions with the 



required vacuum expectation values l|3.10|l below the per- 
cent level, in many cases much smaller. We do not 
find fully breaking global minima for scenarios in the 
range i|3.11[l where the solutions with the required vac- 
uum expectation values (|3.10|) are local minima. Eventu- 
ally, we find examples, where CP conserving parameter 
values with the "wrong" global minimum produce the 
wanted global minimum if a non-vanishing phase S' K is 
introduced. 
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FIG. 2: Same as in Fig. Q but for variation of vs, M H ± , 8' K , respectively. 
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CONCLUSION 



with rational coefficients: 



The stationarity conditions of Higgs potentials at the 
tree-level lead in general to non-linear, multivariate, 
polynomial systems of equations. We have formulated 
the stationarity conditions by means of gauge invari- 
ant functions, thus avoiding to have to deal with elec- 
troweak gauge degrees of freedom. We have introduced 
the Groebner basis approach to solve the systems of 
equations and have applied the method to the NMSSM. 
Within the Groebner basis approach we have easily found 
all stationary solutions even with general CP violating 
parameters. We have demonstrated that large regions of 
the parameter space of the NMSSM are excluded because 
they do not lead to the required global minimum. Finally, 
we note that the method proposed here should also be 
applicable to other potentials which neither have to be 
restricted to two doublets nor even have to be renormal- 
isablc. 
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APPENDIX A: BUCHBERGER ALGORITHM 

In this appendix we want to sketch the construction of 
the Buchberger algorithm which transforms a given set 
of polynomials F into a Groebner basis G. The Groebner 
basis G has exactly the same simultaneous zeros as the 
initial set of polynomials F, but allows better access to 
the actual calculation of these zeros. The general idea is 
to complete the set F by adjoining differences of polyno- 
mials. Before we can present the algorithms themselves 
we have to introduce the two basic ingredients, that is 
Reduction and the S-polynomial. For a more detailed 
discussion we refer the reader to the literature 0, 0, El • 
Here we follow closely |(J. First of all we recall some 
definitions. 

Definition. Polynomial Ring 

A Polynomial Ring K[x±, . . . , x n ] = K[x] is the set of 
all n-variate polynomials with variables Xi,...,x n and 
coefficients in the field K . 

Definition. Generated Ideal 

Let F = {fx, . . . , /„} C K\x\ be finite, F generates an 
ideal defined by 



1(F) 



fiEF 



n e K[x], fi £ F, i = 1, 



In the following we want to consider an explicit exam- 
ple, that is a set F = {fx, / 2 , fs} C Q[x, y] of polynomials 



fx = 3x 2 y + 2xy + y + 9x 2 + 5x - 3, 

f 2 = 2x 3 y -xy-y + 6x 3 - 2x 2 - 3x + 3, 



(Al) 



x 3 y- 



x 2 y 



3x 3 + 2x 2 



The set F generates an ideal 1(F), which is given by the 
set of sums of fx, fi, and fs, where each polynomial is 
multiplied with another arbitrary polynomial from the 
ring Q[a;, y\. The summands of the polynomial are de- 
noted as monomials and each monomial is the product 
of a coefficient and a power product. 

Further, we introduce an ordering (>-) of the monomi- 
als. In the lexicographical ordering (Vi cx ) the monomials 
are ordered with respect to the power of each variable 
subsequently The ring notation Q[x, y] defines y ^i ox x, 
that is for the lexicographical ordering of monomials pow- 
ers of y are considered first, then powers of x. Explicitly, 
this means 2x 2 y 3 >\ Q ^ 5xy 2 because the power of y is 
larger in the first monomial and 2xy 2 >~i ex 5y 2 , because 
both monomials have the same power of y, but the first 
monomial has a larger power of x. The monomials of the 
polynomials IjAlfl from the ring Q[x, y] are ordered with 
respect to lexicographical ordering. In total degree order- 
ing (Vdeg) the monomials arc ordered with respect to the 
sum of powers in each monomial. If two monomials have 
the same sum of powers, they are ordered with respect to 
another ordering, for instance lexicographical. For poly- 
nomials in <Q[x,y] we have x 2 y ^dog 4xy since the sum 
of powers of the left power product is 3 compared to 2 
for the right power product. 

The largest power product with respect to the underly- 
ing ordering (V) of a polynomial / is denoted as the lead- 
ing power product, LP(/), the corresponding coefficient 
as leading coefficient, LC(/). With help of these prepa- 
rations we can define the two essential parts of the Buch- 
berger algorithm, that is Reduction and the S-polynomial. 

Definition. Reduction 

Let f,p G -^[x]. We call f reducible modulo p, if for a 
power product t of f there exists a power product u with 
LP(p) • u — t. Then we say, f reduces to h modulo p, 

tip Coefftcient( f.t) 

where h = j ) ^ J - L 



LC(p) 



■ u - p. 



In the example IjAlf) the polynomial f^ is reducible 
modulo fx, since for example the second monomial of f%, 
that is x 2 y, is a multiple of the LP(/i), and h = fz — 



l/3/i 



x 3 y 



2/3xy- l/3y + 3x 3 



5/3x + l. 



Reduction of a polynomial modulo a set P C ^[x] 
is accordingly defined if there is a p G P such that / 
is reducible modulo p. Further, we say, a polynomial h 
is in reduced form or normal form modulo F, in short 
ixomxf(h, F), if there is no h! such that h reduces to h' 
modulo F. A set P C K[x] is called reduced, if each p G 
P is in reduced form modulo P\{p}. Note that reduction 
is defined with respect to the underlying ordering of the 
monomials, since the leading power product is defined 
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with respect to the ordering. In general, a normal form 
is not unique, neither for a polynomial nor for a set. 

Now we can present an algorithm, to compute a normal 
form Q C K[x] of a finite F C K[x}. 

Algorithm. Normal form 

For a given finite set F C K[x] determine a normal 
form Q C K[x]. 



Q:=F 

while exists p £ Q 

which is reducible modulo Q\{p} do 

Q ~ Q\{p} 

h :— normf (p, Q) 
if h ^ then 
Q:=QL){h} 
return Q 



Algorithm. Buchberger 

For a given finite set F C K[x] determine the Groebner 
basis G C K[x] with 1(F) = 1(G). 



G:=F 

B : = {{9i,92}\gi,92 e G with gi ^ g 2 } 
while B^(3 do 

choose {<?i,<?2} from B 

B:=B\{{g u g 2 }} 

h := spol(gi,g 2 ) 

h' := normf (/i, G) 

if ft,' ^ then 

B:=BU{{ S! / 1 '}| 9 eG} 
G:=GU{/i'} 
return G 



Definition. S-polynomial 

Forgi,g2 £ i^[x] £/ie S-polynomial of gi and gi is defined 



lcm(LP( gi ),LP( g2 )) 
spol(.9i,.92j = \ ffi 



02, 



LP(ffi) 

LC( gi )lcm(LP( gi ),LP( g2 )) 
LC(.g 2 ) LP( 52 ) 

where 1cm denotes the least common multiple. 

In the example I|A1|) we can build the S-polynomial for 
any two polynomials, for instance 

s P ol(/i, h) 4/i-!t/2 =i A- 3 /2 /a 

= 2a; 2 y + 5/2xy + 3/2y + 8x 2 + 3/2x - 9/2. 

Finally we define the Groebner basis. 
Definition. Groebner basis 

G C K[x] is called Groebner Basis, if for all /i,/ 2 £ G 
normf (spol(/i, / 2 ), G) = 0. 

Now we are in a position to present the Buchberger 
algorithm. 



Note, that G generates the same ideal as F, especially, 
both sets have exactly the same simultaneous zeros. It 
can be proven, that the Buchberger algorithm terminates 
as well as that the reduced Groebner basis is unique 0. 
If we apply the Buchberger algorithm to the set IjAlf) with 
subsequent reduction we end up with the reduced Groeb- 
ner basis (with underlying lexicographical ordering) 

gi = y + x 2 — 3/2x — 3, 



fj2 



- 5/2.T 2 - 5/2x. 



The system of equations g\ = g 2 = is equivalent to 
fx = f 2 = fz = 0, but the former allows to directly 
calculate the solutions: Since gi = is univariate it can 
be solved immediately and subsequently g± = for each 
partial solution insertion. 

Despite the correctness of the Buchberger algorithm 
tractability of practical examples requires to improve this 
algorithm. In particular, the number of iterations in the 
algorithm drastically grows with an increasing number of 
polynomials and with higher degrees of the polynomials. 
In this respect much progress has been made with the 
improvement of this original Buchberger algorithm from 
1965 (see [HHH). 
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